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Abstract 



The nearest neighbour level spacing distribution and the A3 statistics of level fluctuations 
associated with very high spin states (/ <; 30) in rare-earth deformed nuclei are analysed by 
means of a cranked shell model. The many particle-many hole configurations created in the 
, rotating Nilsson potential are mixed by the surface-delta two-body residual interaction. The 

£f~^ • levels in the near-yrast region show a Poisson-like level spacing distribution. As the intrinsic 

excitation energy U increases, the level statistics shows a gradual transition from order to chaos, 
£ — , reaching at U ^ 2 MeV the Wigner distribution typical of the Gaussian orthogonal ensemble of 

■ random matrices. This transition is caused by the residual two-body interaction. On the other 

hand, the level spacings between the yrast and the first excited state show a peculiar behaviour, 
displaying a Wigner-like distribution instead of the Poisson-like distribution seen for the other 
near-yrast rotational states. The lowest spacings reflect the properties of the single-particle 
orbits in the mean-field, and are only weakly affected by the residual two-body interaction. 



PACS: 21.10.Re, 21.60.Ev, 23.20.Lv, 24.60.-k, 24.60.Ky, 24.60.Lz 

Keywords: high spin states, level statistics, chaos, cranked Nilsson potential, surface delta interaction 



1 Introduction 

The statistical fluctuations of the energy levels and the transition strengths measured in highly 
excited nuclei with excitation energy above the neutron threshold (several MeV) are well described 
by the random matrix theory |Q |3j. For example, the nearest-neighbour level spacing distribution 
(NND) and the spectral rigidity (or A3 statistics) of the neutron resonance states follow the behaviour 
predicted by the random matrix theory for the Gaussian orthogonal ensemble (GOE) jj], [|. This 
seems to indicate that such excited nuclei, at least over a time scale associated with the observed 
energy interval, are an example of a chaotic quantal system, in the sense that GOE fluctuations 
generally characterize quantum systems which are chaotic in the classical limit |5[ [| . 

The fluctuation properties at lower excitation energy are less well understood, although several 
extensive analysis of low-lying levels as well as of near-yrast high spin levels have been reported 
recently @, |, §, [k], || ||, £§. Although the low lying and low spin levels generally show level 
spacing distributions which are intermediate between chaos (the GOE or Wigner limit) and order 
(the Poisson limit), one observes some systematic behaviour with respect to the mass- number and 
the angular momentum In particular, it is remarkable that the NND in heavy deformed nuclei is 
the closest to the Poisson distribution, not only for the low- lying, low spin levels Q , but also for the 
high spin rotational levels lying near the yrast line O, O, O] . This suggests that both the Poisson 
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and the GOE fluctuations coexist in rotating nuclei and that one should expect a transition from 
order to chaos with increasing intrinsic excitation energy U (the relative excitation energy measured 
from the yrast line at given spin). In the present paper, we examine theoretically the level statistics of 
high spin states in rapidly rotating nuclei as a function of intrinsic excitation energy U. In particular, 
we investigate in detail the level statistics associated with the near-yrast states which may become 
accessible in future experiments. We limit ourselves to the very high spin region with /<;30, where 
static pairing is generally quenched or even vanishes, because our model is not adequate to deal with 
strong pairing correlations. Although this makes it difficult to make a direct comparison with present 
experimental data, there are good reasons to expect that much more experimental information will 
be available in the near future. 

The high spin states near the yrast line in well deformed nuclei form rotational band structures, 
as evidenced by experiments. These rotational band states are usually well described by the cranked 
mean-field models | fl4[ [l5| [l6| in which the collective rotation is represented by uniform rotation 
along the axis of the largest moment of inertia (axis perpendicular to the elongated direction). The 
intrinsic structure of a rotating nucleus is described in terms of the mean-field potential adding the 
cranking term caused by the uniform rotation. Most observed rotational bands are based on intrinsic 
configurations with a few excited quasiparticles (or particles and holes) defined in the cranked mean- 
field Hamiltonian. However, as the intrinsic excitation energy U increases at a given spin, intrinsic 
configurations with many particles and holes (np-nh) will show up and become progressively domi- 
nant. Accordingly, the level density increases significantly, reaching a value around ~ 10 2 /MeV at 
intrinsic excitation energy U ~ 1 MeV above yrast line in rare earth nuclei. One then expects that 
the residual two-body interaction begins to play an important role, mixing the np-nh configurations, 
because the size of its matrix elements (~ 10 keV) is of the same order as the mean level spacing. It 
is also to be noticed that the phenomenon of the rotational damping p7[, which sets in at around 



U ~ 0.8 MeV above yrast 18 , is an important signature of the configuration mixing caused by the 
residual two-body interaction. The fluctuations of the energy levels will be sensitive to the config- 
uration mixing among the np-nh configurations. If the configuration mixing were absent, intrinsic 
excitations would be specified uniquely by the excited particles and holes. In such a situation, the 
level fluctuations may follow the Poisson distribution. On the other hand, once the residual two- 
body interaction is switched on, the np-nh configurations interact with each other. If the residual 
interaction is so strong that many np-nh configurations are admixed with complicated amplitudes, 
one expects that the level fluctuations obey the theory of random matrices. It is therefore important, 
in studying the level fluctuations as a function of intrinsic excitation energy, to take configuration 
mixing explicitly into account. We adopt a shell model approach, making use of a reasonable resid- 



ual two-body interaction on top of a cranked mean field |L9, 20 . Previous work with the cranking 
model has already discussed some general features of the order to chaos transition although it 
used a schematic residual interaction represented by a constant with random sign. We have recently 
shown that the cranked Nilsson model combined with the surface-delta interaction (SDI) |fil| , |22]| can 
reproduce the overall features of rotational damping found in experiment ||^, |j . In the present 
paper we adopt the same model, studying the excited levels lying up to about 2 MeV above yrast 
line. We study in particular detail the states close to yrast, which are likely to be observed in near 
future experiments. 

Statistical analyses of high spin levels in deformed nuclei on the basis of the interacting boson 
model |^{| , the interacting boson fermion model p3] , and the particle-rotor model p7l have also been 
reported. These models, however, take into account only limited degrees of freedom (sd collective 
bosons or high-j nucleons) of the intrinsic excitations in deformed rotating nuclei. 

2 Formulation 

2.1 The model 

We start with the cranked Nilsson single-particle Hamiltonian 

hcrank — ^Nilsson ^jx (1) 
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in order to define the single-particle basis in a rotating deformed nucleus. Here the quadrupole and 
hexadecapole deformations are considered. We do not include the static pairing potential in the 
mean-field. This may be justified for the high spin region (I 30) which we are mostly concerned 
with, since the pairing gap is usually reduced, or even vanishes, due to the rotational perturbation 
(Mottelson-Valatin effect) ^8|, |2^, |30| . The eigen-solutions of the cranked Nilsson single-particle 
Hamitonian define an adiabatic basis as a function of the rotational frequency to. However, since the 
adiabatic orbits sometimes accompany avoided crossings between orbits which cause abrupt change 
of the basis wave functions against small change inw, we instead use a diabatic single-particle basis, 
which is constructed by removing small interactions causing the repulsions at the avoided crossings. 
Putting N neutrons and Z protons in the diabatic single-particle basis, shell model many-body 
configurations (labeled by /i) are generated: 

im(/)>= n a ii->- (2) 

occupied i in fj. 

In Eq.(2) aj denotes the nucleon creation operator for an occupied diabatic single-particle orbit i, 
which is defined at an average rotational frequency uji corresponding to the given angular momentum 
I. We include all the single-particle orbits within an interval of 3.0 MeV below and above the Fermi 
surface. The shell model basis includes the configuration in which the single-particle orbits 

up to the Fermi surface are fully occupied, as well as all possible np-nh configurations with respect 
to the fully occupied one. 

The energy of a shell model configuration is given, following the standard cranked Nilsson- 

Strutinsky prescription, by 

E„(I) = E^ tls (I) - E smooth {I) + E RLD {I) (3) 

where E^ lls (I) = E'^lo) + loJ X: ^(lo) with the angular momentum constraint J X41 (lo) = I on the 
rotational frequency u>. Here E'^(uj) = J2ien e i( w ) an< ^ = Sieui^.i^) are the total routhian 

and the expectation value of the angular momentum J x of the shell model basis /x, respectively. Since 
we use the diabatic single-particle basis which depends only weakly on the rotational frequency, the 
energy expression can be accurately approximated locally by 

E™°(I) = + UI I + (/ ' Ja: ^ j))2 (4) 

2 J ^ 

(2) 

referring to the average rotational frequency uji. Here Jjj, is the dynamical moment of inertia of 
the configuration. The deviation | J XiAI (w/) — I\ in the angular momentum expectation value is less 
than 5 at spin 1 — 50 for most configurations in the present calculation. Although the Strutinsky 
smoothed energy E smooth (I) and the rotating liquid drop energy E RLD (I) correct the absolute 
excitation energy, they do not affect the level statistics discussed in the present paper. 

We then introduce a two-body force, mixing the shell-model configurations. We adopt the surface 
delta interaction (SDI)Q 

v(l, 2) an ^ = -4nV ]T Yl M (d t ,i(l>t,i)YLM{e t ,2(t>t,2) (5) 

LM 

where {0t4>t) is the angle variable in the stretched coordinates. The strength parameter Vo includes 
the radial matrix elements and we use the strength Vo = 27. 5/ A MeV given by Ref.[^2|, which is the 
same value used for the study of rotational damping in 168 Yb |2(], |23| . The shell model Hamiltonian 
is given by 

#(/W = -e^OO-W + v(iu> (6) 

where T^(/) MM ' denotes the matrix elements of the residual two-body interaction of SDL The Hamil- 
tonian is diagonalized to obtain energy eigenstates 

\a{I))=J2 X "(I)W)) (7) 
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K>(U6l Dy 161,l(j2jj 
162,163Qy 163,164jj q 
164,165Qy 165,166jj q 
162,163g r 163,164rp m 
164,165g r 165,166rp m 
166,167g r 167,168rp m 
166,167y]-, 167,168l u 

i68,i69Yb i69,i70j jU 

170,17lYb 171 > 172 Lu 
172.173Y13 173 . 174 Lu 
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0.248 
0.261 
0.267 
0.245 
0.258 
0.267 
0.246 
0.255 
0.265 
0.269 



-0.016 
-0.007 
0.003 
-0.009 
0.001 
0.012 
0.004 
0.014 
0.025 
0.036 



Table 1: The quadrupole and hexadecapole deformation parameters 62 and £4 used in the present 
calculations 



which are admixtures of the basis configurations {\p(I))} as well as their energy levels {E a (I)}. The 
diagonalization is done separately for each I n , truncating the basis by including the lowest 1000 
basis states. The resulting lowest 300 states (covering the region up to U ~ 2.4 MeV) are rather 
stable against the truncation of the basis. For further details of the model, we refer to Ref. p(J. 

In the present paper, we focus on the rare-earth nuclei, and in particular we consider 40 nuclei 
in the A = 160 — 174 region, listed in Table 1, for which deformed prolate shape stable up to very 
high spins is suggested by potential energy surface calculations [|l[ We adopt the equilibrium 
deformation parameters taken from Ref. |33]| , as given in Table 1, which are similar to those calculated 
in Ref. |Tl| . In order to make a statistically meaningful analysis, we collect the spacings taken from a 
certain spin interval in all the 40 nuclei, and we will not discuss the dependence on individual nuclei, 
spins, and parities. 

In the following, we use the parity and the signature quantum number (tt, a) to classify the 
energy levels of the total system. The signature a is related to the total spin / through the relation 
I = Io + a with a = 0, 1 for even- A system, and a = ±1/2 for odd- A We sometimes use the even 
integer spin Iq and the signature a in place of the "true" spin / when we specify spin intervals. 

2.2 Level statistics 

In order to perform the statistical analysis of the energy level fluctuations, one must take into account 
the fact that the level density and hence the level spacing are strongly dependent on the intrinsic 
excitation energy U. In this situation, it is necessary to separate local level fluctuation from the 
overall excitation energy dependence of the level spacings. For that purpose, we adopt the unfolding 
procedure | in a particular form which follows Shriner et.al. [^|. The unfolding procedure 
measures the local level fluctuations with respect to a smooth average level density. We assume that 
the average level density is represented by the constant temperature formula |35[] 



P(£)^exp(^) (8) 

for each spectrum at a given I 71 . To determine the parameters in the formula, we make a fit to the 
staircase function which represents the cumulative number of levels below energy E, 

/E 
p{E')dE' = Y,Q{E~E a ) , (9) 
-00 a 

with a smooth function corresponding to the average level density p(E) , 

'E-Ea 



N(E) = / p(E')dE' + N a = cxp ^ - ] - 1 + No , ( LI) i 



4 



by minimizing the quantity 



G(T,E Q ,N ) = (N(E)~N{E)) 2 dE (11) 

J E m i n 

with respect to the parameters T,Eq and Nq for each spectrum at a given I 71 ". Here the energy 
boundaries E m i n and E max are the energies of the lowest and the 300-th levels in each spectrum. 
When we discuss the level statistics for the lowest 20 levels, however, we obtain a better fit including 
only the lowest 30 levels. The unfolded spectra {x a ; a = 1, 2, ...} are then derived for each I" by the 
transformation 

x a = N(E a ). (12) 

The unfolded spectra have a constant average density p x (x) = 1 provided that the constant temper- 
ature formula fits well the average level density. 

In order to analyze the level fluctuations, we calculate the nearest neighbour level spacing distri- 
bution (NND) which is also often used in the experimental analysis. We calculate the distribution 
P(s) for the unfolded spectra, where s = x a +i — x a is the spacing between the neighbouring levels 
with same T" . By the unfolding procedure, the spacings are normalized as (s) = 1. The distribution 
is represented as a histogram. 

The NND is calculated for various ensembles of level spacings which are taken from different 
intervals in excitation energy. The obtained distribution is fitted with the Brody distribution jl] 

l+tu 



P w (s) = (l + w)as w exp(-as L+w ), a= \T\^——JJ , (13) 

parametrized by the Brody parameter w. This family of distributions is convenient because the 
Brody parameter w = 1 produces the Wigner distribution while the value w = corresponds to the 
Poisson distribution. (Note that the theory of GOE random matrices leads to w — 0.953 [Q , which 
is not distinguishable from the Wigner limit in the present analysis). The value of w is determined 
minimizing the quantity 

*w = E(^P) 2 . 

where P(i) = N(i)/N is the probability in the i-th bin [si,Si + As] of the calculated NND (N(i) 
being the number of spacings in bin i out of the total spacings of N). P w (i) — f*' s P w (s')ds' 
is the corresponding probability in the Brody distribution. The statistical error is estimated as 
cr(i) — sj N(i)/N for N(i) > and a(i) = 1.15/N for N(i) — by assuming the multinomial 
distributions. 

We also calculate the ensemble average of the A3 statistics J3(| [l| 
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rx+L 

A 3 (L) = ( j- min / [N x (x') - Ax' - B] 2 dx' ) (15) 



or the spectral rigidity. Here N x (x) — *^2 a 0(x — x a ) is the staircase function for the unfolded 
spectra and the average (...) is calculated over spectra in a given ensemble and intervals [x,x + L], 
[x + L/2, x + 3L/2], ... in a spectrum ||. For the Poisson distribution of levels, 



while for the GOE distribution. 



A 3 ,Poi S so„(£) = L/15 , (16) 



A 3 ,GOE(i) « -4(lnL - 0.0687) . (17) 
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Figure 1: The Brody parameter extracted from the NND for energy bins containing the first to 
5-th, 6-th to fO-th, 11-th to 20-th, 21-st to 30th, and 31-st to 40-th levels, ... 291-st to 300-th 
of each spectrum. The result is plotted as a function of the intrinsic excitation energy covered 
by the bins. The solid, dotted, and dot-dashed lines correspond to different spin intervals la = 
32 - 50, 20 - 30, 52 - 60 respectively. 



3 Results and Discussion 
3.1 Order to chaos transition 

We first discuss how the level statistics depends on the intrinsic excitation energy U, aiming at 
extracting the overall dependence on U in a wide interval ranging from U = (at yrast line) to 
U ~ 2 MeV. For that purpose, we calculate the NND and A3 for the lowest 300 levels in each 
spectrum, grouping them in bins of levels. The intrinsic excitation energy of the binned levels 
approximately covers the region up to U ~ 2.4 MeV. 

The calculated Brody parameter for the NND for the spin interval I = 32 — 50 is depicted 
in Fig|l]. The Brody parameter increases monotonically with increasing intrinsic excitation energy 
U. The NND for the lowest bin (first to 5-th levels at each /*'), has the Brody parameter w — 
0.301 ± 0.012. The corresponding NND shown in Fig^(a) is much closer to the Poisson than to the 
Wigncr distribution, although one can also notice a small deviation from the Poisson distribution. 
For the levels from 10-th to 40-th, the Brody parameter is about 0.5, which is midway between the 
Poisson and the Wigner distributions, as can also be seen from the NND plotted in Fig^(b). As the 
intrinsic excitation energy increases further, the NND approaches the Wigner limit; the bin including 
the levels from the 291-st to 300-th has w — 0.888 ±0.012, which is close to the Wigner limit w — 1.0 
or GOE limit 0.953 (See also the NND shown in FigJ|(c)). These results indicate that the transition 
from order (Poisson fluctuation of the levels) to chaos (Wigner and GOE fluctuations) takes place 
gradually increasing the intrinsic excitation energy, until the chaotic limit is nearly achieved at 
around U ~ 2 MeV above yrast line. This dependence on excitation energy confirms the results of 
a previous analysis p3| performed with the same model, but without transforming the energy from 
the rotating frame to the laboratory frame (the last two terms in Eq.(Q) were neglected). 

It is interesting to note the implications of the results from the A3-statistics (Fig||). The GOE 
limit is obtained only for L values up to some value L m ax, and it is found that L max increases with 
increasing excitation energy. For the bin at the highest studied excitation energy (#251 — #300, 
U ~ 2.3 MeV), we find that L max ~ 6. This implies that an energy eigenstate in this interval follows 
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the GOE correlation only with approximately the ten closest lying states. Thus, the GOE-behaviour 
seen in the NND for the energy levels in this interval (FigJ|(c)) indicates a chaotic behaviour of 
local nature. The NND and A3 carry different types of information for short-range and long-range 
correlations, as discussed in (37). Since the spreading width of np-nh shell model basis states is 
finite, L max could be related to pT^ (p being the level density) jl{| [37] although an estimate based 
on L max ~ 2.hpT ^ which is found in a random matrix model |37j gives a much larger value than the 
calculated L max ~ 6. Non generic behaviours of A3 have also been discussed in connection with the 
shortest periodic orbits p8| and the Lyapunov exponent |39j in semiclassical analysis, whose relation 
to the present model is not clear yet. 

In Fig.|l|, we also show the Brody parameters extracted from ensembles of binned levels taken from 
lower and higher spin intervals with Iq = 20 — 30 and Iq — 52 — 60. No significant spin dependence 
is observed. This is in contrast with the interacting boson fermion model p6j and the particle-rotor 
model |^7j], which predict a spin dependence caused by the alignment of the high-j orbitals. Note 
that, besides the high-j orbitals, we include all the other single-particle orbits near the Fermi surface, 
which do not necessarily align in the considered spin interval. 

It is interesting to compare our results with the previous theoretical analysis by Aberg |l9| , who 
used essentially the same model except for the matrix elements of the two-body interaction, which 
were schematically approximated by a constant with random sign. When the mean square root 
value of the matrix elements was 15 keV, the A3 statistics reached the GOE limit in the excitation 
energy range C7 = 1.5 — 2.0 MeV in 168 Yb , a value lower than in the present model. The difference 
can be traced back to the statistical properties of the two-body matrix elements. We find that the 
statistical distribution of the off-diagonal matrix elements of the SDI force follows a distribution 
strongly peaked at zero matrix element compared with the Gaussian distribution [ 
a selectivity in the two-body matrix elements related to the intrinsic nature of the SDI and of the 
cranked Nilsson single-particle orbits. Because of this selectivity, the onset of chaos in our model 
takes place at higher excitation energy, although the average mean square root of the off-diagonal 
matrix elements is about 19 keV in the present calculation. 



20. 23), indicating 



3.2 Level statistics in the near-yrast region 

The bin including the lowest 5 levels for each I* covers the interval up to U ~ 0.7 MeV above 
the yrast line. The calculated levels in this energy region mostly form rotational band structures 
connected by strong stretched E2 transitions [M. These levels are probably those which will be 
resolved in experiments in the near future, while it will be much harder to resolve excited levels lying 
in the region of rotational damping (U <; 0.8 MeV). In fact, up to around 10-20 rotational bands are 
observed in a few rare-earth nuclei j|0[ ^] (although only a few rotational bands are identified at 
the highest spin / ~ 40). In this subsection, we discuss in detail the level statistics associated with 
these near-yrast states. 

For this purpose, we introduce a strict ordering of the spacings according to excitation energy 
above yrast. The strict ordering N encompasses four spectra having different parity tt and sig- 
nature a for a given reference spin Iq (even integer), that is those spectra with I v = Iq,(Iq + 
1) ± for even- A systems, and those with i" w = (Iq ± l/2) ± for odd- A systems. More precisely, 
we first define a reference energy E re f(I) by an envelope of the yrast levels, i.e, the E re f(I) — 
min (Ei owes t(I), - E '°"' e3t ( J + 1 )+- E '°'" e3f ( J ~ 1 ) j m orc [ er to compare the four spectra. We then assign 

the label N to the levels in the four spectra, counting from the lowest according to the excitation 
energy E(I) — E re f(I) measured from the reference. By this definition the N = 1 level represents 
the "strictly yrast" level in the sense that it refers to only one among the four lowest levels defined 
separately for the four parity and signature quantum numbers, while the other three are treated 
as excited levels (N > 1) with respect to the strict yrast level. Collecting the spacings from the 
iV-th rotational band (they are the spacings between the iV-th level and the next excited level with 
the same 1^) from the spin range Jo = 32 — 50 in the 40 rare-earth nuclei, we calculate the NND 
and extract the Brody parameter for each N. We also made the same analysis for spin intervals, 
Iq = 20 — 30 and Iq = 52 — 60 in order to study a possible spin-dependence, although the present 
model may not be very realistic for the lower spin interval (Iq = 20 — 30) because of the problem of 
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Figure 4: The Brody parameter extracted from the NND associated with the lowest 15 near-yrast 
states. See text for the definition of the strict ordering N of the states. Different symbols represent 
spin intervals Iq = 32 — 50, Iq = 20 — 30, and Iq — 52 — 60. 

the pairing correlation. 

The extracted Brody parameter is plotted in Fig. | as a function of the strict level order N, 
for the lowest 15 levels. The spacings analysed for Fig. ^ belong mostly within the first energy 
bin adopted in the previous subsection, which included the lowest 5 levels at each I w . The Brody 
parameters plotted in Fig. ^ in average agree with the value w ~ 0.3 for the lowest bin in Fig.|l|. 
It is seen in Figj^ that the Brody parameter gradually decreases as the levels become closer to the 
yrast line, except for N—l. For the lowest few levels, the Brody parameter is about 0.1-0.2, which is 
close to the Poisson limit (The corresponding NND's are shown in Fig.||(b,c) for the third and tenth 
lowest levels N — 3, and 10). This indicates that the excitation energy dependence shown in Fig.|l| 
continues down to the lowest few states near the yrast line. However, a remarkable deviation from 
the overall excitation energy dependence is clearly noticed for the lowest point at N = 1, i.e., for the 
spacings between the yrast band and the next excited band with same I v , for which w ~ 0.4 — 0.7. 
The NND for N = 1 is shown in Fig|(a). 

It is also seen that the spin dependence is not strong while at N = 1 the Brody parameter at 
lower spins shows a more significant deviation from the Poisson limit, becoming close to the Wigner 
limit. 

3.3 Level spacing statistics at yrast 

In order to study the origin of the deviation of the first spacings from the Poisson distribution, we 
perform calculations neglecting the residual two-body interaction in which all the states become pure 
many-particle many- hole mean-field configurations. The NND and the extracted Brody parameter 
are compared in Figs.|| and |^ with those obtained by inclusion of the residual interaction. 

The Brody parameter of the lowest spacings (N = 1) is essentially unaffected by the inclusion 
of the residual interaction. This indicates that the deviation from the Poisson limit for N = 1 does 
not originate from the residual interaction, but from the mean field. On the other hand, Figures ^| 
and ^ also show that for the higher spacings above N ~ 3, the Brody parameter of the pure mean 
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I =32-50 




Figure 5: The NND associated with the lowest, third, and the tenth (strict order 7V=1,3,10) levels 
within the spin interval Iq = 32 — 50 for (a), (b), and (c), respectively. (d,e,f) the same as (a,b,c) 
except that the residual SDI interaction is neglected. 
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Figure 6: The Brody parameter extracted from the NND associated with the lowest 15 near-yrast 
states in the strict order obtained from the mean-field calculation without the residual interaction 
for different spin intervals Iq = 20 — 30, Iq — 32 — 50 and Iq = 52 — 60 (points joined by dashed 
curves). It is compared with the result with the residual interaction (ones joined by solid curves. 
See also Fig.|). 
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field calculation converges to the Poisson limit w = except for Iq — 20 — 30 [] while the Brody 
parameter and the NND calculated with the residual two-body interaction deviates from the Poisson 
limit. This explicitly indicates that , contrary to the case of N = 1, the deviation from the Poisson 
limit at N <; 2 arises from the residual two-body interaction. 

We look for the origin of the special feature of the very yrast N = 1 spacings in connection with 
the single-particle level structure in the cranked Nilsson mean-field. To this end, we first remark that 
for the levels near the yrast line, the configurations are only weakly mixed by the residual two-body 
force and most of them have essentially independent-particle configurations. In particular, the excited 
states near the yrast lines often have lplh configuration with respect to the yrast configuration. 
This means that the relatively large value of the Brody parameter extracted for the yrast levels may 
not be directly related to the mixture caused by the residual interaction. Furthermore, when the 
spin is not very large, the angular momentum alignments of intrinsic excitation is relatively small 
compared to the level spacings between the states in the yrast band and the next excited states 
with the same quantum number; it is found that the last term in Eq.(|J) representing the alignment 
effect on the energy is at least a factor ^ 2 smaller than the average level spacing D ~ 350 keV 
associated with the N = 1 yrast rotational band. Under these conditions, the relative energy of the 
first excited states having the same J 7r measured from the yrast states can be approximated by the 
lplh excitation energy in the single-particle routhian spectrum. Namely, E next (I n ) — E yraa t{_I' K ) ~ 
e' p alT {u>i) — e' h a7V (u>i) where e' p h are the single-particle routhian of the involved particle and hole. 
The particle and hole orbits necessarily have the same quantum numbers (signature a and parity 
7r). On the other hand, spacings between excited states (yrare levels) do not keep such relation to 
the single particle spacings. This is illustrated in Fig. [7] which shows examples of the main mean 
field components of the lowest excitation within each parity-signature set of states. It is clearly seen 
how the excitation from the yrast state in this typical case will proceed by changing the orbit of one 
particle, keeping its (it, a) (as shown in the left panel of Fig.^). For yrare levels, excitations relative 
to the Fermi surface are already present, and the lowest excitation starting from an yrare state most 
often will proceed by a 2p-2h excitation which connect orbitals with different (tt, a), as shown in the 
right panel. 

In the previous section, we find that the the special feature seen for the yrast bands becomes 
more prominent as the spin decreases. This feature is consistent with the above interpretation. In 
addition, we find that there exists an odd-even effect for N ^ 10 for the low spins Iq = 20 — 30; 
the level statistics for odd-odd nuclei is very close to the Poisson distribution while even-even or 
odd- A nuclei shows significant deviation from the Poisson distribution for N — 1 . This also indicates 
that the level spacings associated with the very yrast states reflects the single-particle level spacings 
since the odd-even effects can arise from the fact that many cranked-Nilsson routhian orbits retains 
two-fold degeneracy (signature splitting is small) at low rotational frequency. 

These considerations lead us to investigate the spacing distribution of the single-particle levels in 
the cranked Nilsson potential. Figure || shows distribution of the level spacings e' pan (u>i) — e' h an (wi) 
between the hole and particle orbitals having the same parity and signature quantum number which 
correspond to the lowest lplh excitations for all 40 nuclei and rotational frequencies corresponding to 
spin I — 30, 32, ...50. We have not applied the unfolding procedure since the relevant single-particle 
orbits lie only in a limited region around the neutron and proton Fermi surfaces (N = 94 — 105, Z = 
66 — 71) of the cranked Nilsson spectrum. It is noticed in Fig.|| that the spacing distribution is 
concentrated around the average spacing ((D) ~ 400 keV) and there are few spacings smaller than 
200 keV, indicating that degeneracy among orbits with the same quantum numbers happens only 
rarely. This in fact comes from the nature of the cranked Nilsson spectrum (which is believed to 
be valid also for other models such as Woods-Saxon potential). One of the relevant properties is 
that a large part of the deformed mean field is of harmonic oscillator type, and that the quantum 
spectrum of the oscillator shows strong level repulsion while the corresponding classical motion is 
integrable. With deformation not very large (e Si 0.3) the single-particle orbits around the Fermi 

For the rotational frequencies corresponding to the spin Io = 20 — 30, many of the cranked-Nilsson routhian orbits 
show very little signature splitting. This causes frequent near degeneracy in the rap-nh configurations, leading to 
the enhancement of the NND for small spacings and producing negative values of the Brody parameter for N J> 3 in 
I = 20 - 30 case in Fig.|. 
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Figure 7: Main mean-field configurations in the lowest states with (n, a) = (—,1) and (+,0) in 
168 Yb at w = 0.319 corresponding to I = 30,31. The blobs represent the occupied orbitals in 
the main configuration of the lowest state for each (tt, a). At this spin, the yrast state belongs to 
(— ,1). The arrows represent excitations involved in the second lowest state in each (n,a). The 
solid, dotted, dashed and dot-dashed lines represent the cranked Nilsson single-particle orbits with 
(tt, a) = (+, 1/2), (+, -1/2), (-, 1/2), and (-, -1/2) respectively. 



surface belong to a single major oscillator shell, provided that the parity and the kind of particle 
is fixed. Taking the neutron spectrum as an example, the negative parity orbits are dominated by 
those with the total oscillator quanta iV osc = 5. Because of the mean- field deformation, the orbits 
having different 713 (oscillator quanta along the deformation axis) are then splitted in energy, and this 
makes degeneracy among orbits with different 113 asymptotic number rare. Furthermore, the I 2 and 
Is terms of the mean field cause splittings among the orbits having the same 72,3. The positive parity 
neutron orbits near the Fermi surface are £13/2 orbits, and because of the deformation splitting, 
the ?i3/2 orbits with fixed signature are placed with finite intervals at any rotational frequency. 
Therefore degeneracy among the 113/2 orbits never happens. An additional mechanism arises from 
the fact that the nuclear mean-field favours an equilibrium deformation at which the shell energy 
lowers, implying that degeneracy of single-particle orbits at the Fermi surface is unfavoured. All 
these mechanisms prefer to the Wigner-like distribution in the single-particle spacing at the Fermi 
surface. Consequently, there exist only few cases of small spacings as seen in Fig.||, especially for the 
low spin region Iq < 30. At higher spins (i.e., at high rotational frequency), some specific orbits with 
large rotational alignment, e.g, proton orbits stemming from hg/2 and I13/2 intrudes in the Fermi 
surface region around I ~ 40 — 50 and cross sharply with other orbits. Small spacings associated 
with these highly aligned orbits are present in the distributions shown in FigJ|(b) for Iq = 32 — 50 
(and slightly also for Iq = 52 — 60), but this does not enhance very much the probability of small 
spacings and keep the distribution Wigner-like. 

Consequently, the intrinsic nature of the cranked single-particle spectrum affects specifically the 
level spacing distribution associated with the yrast band at N = 1. Figure |6| indicates that this 
remains to some extent even for the very high spin region with / = 30 — 60. By the same token, the 
present analysis suggests that the singular behaviour of the lowest spacing could be stronger at lower 
spins. It should be remarked however that the present model is not very accurate for describing 
the near-yrast rotational bands at lower spins since the pairing correlation which is important at 
low spins is not well taken into account. Thus, the results of the present analysis cannot readily be 
compared to experiments for the lower spins (/ 30) . 
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Figure 8: (a) The distribution of the spacings of the cranked Nilsson single-particle orbits for 
spin interval Iq = 20 — 60. The sampling is described in the text. The dotted line represents the 
distribution which is obtained when the adiabatic basis is adopted. There is no significant difference 
between the diabatic and adiabatic basis, (b) The same as (a), but the histogram bins are defined 
with the spacing itself instead of the normalized spacing, and also the spin interval is subdivided 
into I Q = 20 - 30, J = 32 - 50 and I = 52 - 60. 



We have here related the angular momentum dependence of the spacings of the single-particle 
spectrum in the mean field, displayed in Fig.|^(b), to the behavior of specific important orbits in 
the cranked potential. One may ask whether a more consistent explanation may be achieved from 
the general properties of the phase space associated with the classical single-particle motion in 
the rotating potential. So far, the questions of chaotic and regular motion in rotating deformed 
potentials have only been carried out for billiards in two dimensions |4^, It is found that 
rotation certainly may affect the phase space generally. However, especially for the high kinetic 
energies of nucleon states around the Fermi surface, the predicted effects are small. 

The absolute nature of the yrast band relative to yrare bands of the other parity-signature 
configurations is worth emphasizing. We have discussed it from our model, especially by means 
of FigJ^. However, it may actually be less specific to the actual model, since it could result from 
a general homo-lumo gap in a quantal system, but now for states within an interval of angular 
momenta. The favoured yrast configuration should then be able to determine the detailed shape and 
other properties of the nucleus, and the yrare levels then have to adjust to this. 



3.4 Analysis without unfolding procedure 

As described in Sect.^J, we apply the unfolding procedure in order to separate the overall excitation 
energy dependence and the local fluctuations in the level spacings. This procedure, however, cannot 
be used for the analysis of the present experimental data in the high spin region. In fact, the number 
of identified levels at fixed I™ is far below 10 at high spins in the experiments performed so far. 



The experimental analysis in Refs. |jll|, 12, [L3| does not apply such an unfolding procedure that is 
described in SectJ|. In order to facilitate a direct comparison between the theoretical results and 
experiments, we propose in this subsection another way of analysis which does not use the unfolding 
procedure, but still takes into account the excitation energy dependence of the level density in an 
approximate way. The procedure is also applicable to the analysis of experimental data. 

Although the level density increases exponentially with increasing intrinsic excitation energy U, 
it may be assumed that the level density at given U is rather independent of spin and parity and 
nuclear species as far as the high spin states in the same mass region are concerned. In fact, as 
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Figure 9: The Brody parameter extracted from the NND obtained by using the simple normalization 
without the unfolding procedure for the spin interval Iq — 32 — 50 as a function of the strict level 
ordering N (the symbols connected with the solid line), compared with that obtained with the 
unfolding (dotted line). The inset shows the NND for N = 1 obtained without unfolding. 



discussed in Ref. [Q , the level density at fixed I* can be accounted by the level density of intrinsic 
configurations in the cranking model if the spin is sufficiently high (e.g. above I i5 12 for U < 3 MeV). 
In this limit, the level density may be approximated by the Fermi gas formula p4[ ] as a function of a 
single variable all where a is the level density parameter related to the single-particle level density 
at the Fermi surface of the cranked mean-field. 

Keeping this in mind, let us consider an ensemble of level spacings which is specified by the 
strict level ordering N as introduced in the previous subsections. The level spacings within this 
ensemble are expected to have a common average value since the level ordering N and hence the 
excitation energy is taken to be the same. It then may be reasonable to define, without using the 
unfolding procedure, the normalized spacing s = D/(D) by simply dividing the spacing D by the 
average spacing (D) calculated for this ensemble specified by N . We show in Fig.|| the NND and the 
extracted Brody parameter calculated in this way. It is seen that there exits small but systematic 
difference for N 5s 5 between the results calculated with and without the unfolding procedure. The 
origin of the difference can be understood by noting that the average level spacings for the lowest 
states are 386, 271, 214, .. 163 keV for N = 1,2,3, ..5, which are not very small compared with 
the temperature parameter T ~ 350 keV in the fitted level density. In other words, the smooth 
level density p(E) varies significantly in the energy interval of the single spacing, especially if the 
lowest few N's are concerned. This causes a difference in the profile of the NND depending on 
whether we adopt the unfolding procedure or not. However, it should be stressed that, in spite of 
the difference depending on the way of analysis, the Wigner-like property associated with the yrast 
spacings (N = 1) is present in both analysis. It becomes even more significant with the simple way 
of analysis without the unfolding procedure. 

Except for the lowest several spacings in the strict ordering, the NND obtained by means of 
the simple normalization agrees with those obtained with the unfolding procedure. As another 
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Figure 10: The Brody parameter extracted from the NND associated with the n-th level at each I* 
for spin interval Iq = 32 — 50 (See text). Here the unfolding procedure is not applied. The result is 
compared with the one with the unfolding procedure which is calculated for the bins of levels (same 
as the solid line in Figj^). 



illustration, we consider an ensemble of level spacings specified by the level ordering n defined in 
each spectrum for each /*' (note the difference between n and the strict ordering N) , and calculate the 
NND and the Brody parameter for spin interval I — 32 — 50 by means of the simple normalization 
procedure described above. In this case, the average spacing (D) is calculated for each n. The 
result is compared in Fig.|lO] with the Brody parameter (Figj^ in subsect.3.1) analysed for level bins 
n = 1 — 5,5 — 10, 11 — 20, ... by using the unfolding procedure. Both ways of analysis agree very well 
with each other, leading to the same conclusion about the overall excitation energy dependence of 
the NND. The NND's for n = 1 and 2 are plotted in FigQa,b). 



3.5 Relation to experimental analysis 

The experimental NND obtained by Shriner et al. Q for the low-lying low-spin states (/ S5 5h) in 
rare-earth nuclei displays a Brody parameter around 0.3. A more recent analysis which includes the 
observed rotational states at relatively high spins (most of the analysed levels have / 30) reports a 
NND which is close to the Poisson distribution Jll], [l^, |l3). Our theoretical calculations for higher 
spins / iS 30 also favours a Poisson-like NND for the levels near the yrast line. In the following we 
try to perform our theoretical analysis in a way similar to the procedure adopted by Garrett et al. 
]lT| . [l2] , |l3| . One should however remark that a comparison between our results and the experimental 
findings can only be indicative, because they refer to different spin regions. 

In accordance with Ref. (jll]], we consider here an ensemble of the spacings associated with the 
lowest and second lowest states at each J 7r (n = 1 and 2). However we deal with the spin interval 
Io = 32 — 50, where the pairing effects are expected to be weak. We remark that the pairing effects 
are removed to some extent from the experimental analysis by excl uding the lowest (0,+) spacings 
fTl| . We adopt the normalization scheme introduced in Subsect. 3.4 which do not use the unfolding 
procedure since the experimental analysis |p"Tf adopt the similar normalization. The obtained NND's, 
shown in Fig.[ll|(a,b), are close to the Poisson distribution, having Brody parameter w ~ 0.25 in 
agreement with the Poisson-like NND seen already for the near-yrast states. It also shares some 
common features with the experimental analysis fll|| : A deviation from the Poisson distribution is 
seen for small spacings s Ss 0.2 and is most significant for the smallest spacing with s 0.1 or D ^ 25 
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Figure 11: The Brody parameter extracted from the NND associated with the lowest and second 
lowest states (n = 1, 2) at every P r , for (a) and (b), respectively. Here the unfolding procedure is not 
applied, and the spacings are normalized to the average spacing defined in each ensemble, (c) and 
(d) , the same as (a) and (b) except that the residual SDI interaction is neglected and pure mean- field 
many-particle many-hole configurations are considered. 



keV. 

In Fig.0(c,d), we also calculated the NND in the same way except that the residual two-body 
interaction is neglected. Comparing Fig.[il|(a,b) and Fig.|li"|(c,d), it is indicated that the deviation 
from the Poisson limit at small spacings, seen in Fig.[ll](a,b), is mostly caused by the residual two- 
body interaction. We remark, however, that the Wigner-like distribution associated with the very 
lowest spacing N = 1 discussed in Subsect. 3.3 should be present, but is not visible neither in 
Fig.[ll](a) nor in (c). This is because the ensemble with n = 1 contains also the other spacings 
with N = 2,3, .., which have lower average spacing and mask the Wigner-like distribution associated 
with N = 1 spacings. This suggests that, in order to find the Wigner-like distribution caused by 
the mean-field effect in the experimental analysis, one should make an analysis by subdividing the 
ensemble of the spacings with respect to the strict level order N. 



4 Conclusions 

We analysed the level statistics of the high spin states with I <; 30 in rare-earth deformed nuclei as 
a function of the intrinsic excitation energy of the rotating nuclei. We used a shell model approach 
which describes np-nh excitations in the cranked Nilsson potential interacting through the surface- 
delta residual interaction. We put emphasis on the analysis of the near-yrast levels which may be 
accessible by discrete 7-ray spectroscopies. 

The nearest neighbour level spacing distribution (NND) and the A3 statistics indicate that the 
level fluctuations in the near-yrast region follow a distribution close to the Poisson limits with the 
extracted Brody parameter w — 0.2 — 0.3. This value of the Brody parameter implies significant 
deviation from the Poisson distribution for spacing smaller than about s^0.3. The experimental 
analysis of the NND for low-spin states [|| and the one including high spin rotational states Jy], [l2|, [H| 
indicates a Poisson-like distribution in the near-yrast states. The present analysis suggests that this 
behaviour extends up to very high spins I ~ 50 — 60. The level statistics approaches the GOE limit 
as the intrinsic excitation energy U increases, but this process proceeds very gradually and the chaos 
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limit is nearly attained only with £7^2 MeV. This transition is caused by the residual two-body 
interaction. 

An interesting aspect of the NND emerges when we focus on the lowest levels near the yrast line. 
The level spacings between the yrast rotational band and the next excited band with the same spin 
and signature favour a Wigner-like NND, rather than obey the Poisson-like distribution associated 
with the other near-yrast levels. The distinguishable property of the NND associated with the 
yrast rotational bands arises from the mean-field property of the rotating nuclei while the deviation 
from the Poisson limit seen for the other spacings among yrare rotational bands is caused by the 
residual two-body interaction. Since the lowest few levels near yrast have dominantly one particle 
excitations, the spacing between the yrast rotational level and the next excited level thus reflects 
the single-particle routhian spectrum in the rotating mean-field at equilibrium normal deformation, 
which shows a Wigner-like distribution for the spacings between orbits with the same parity and 
signature around the Fermi surface. 
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